Disclination Dipoles as the Basic Structural Elements of Dielectric Glasses 
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We show that the experimentally observed behavior of thermal conductivity of dielectric glasses over 
a wide temperature range can be explained by a combination of two scattering processes. The first 
one comes from the phonon scattering due to biaxial dipoles of wedge disclinations while the second 
one is the Rayleigh type scattering. The results obtained support the cluster picture suggested 
earlier for glassy materials. 



Physics of glassy systems has been attracting a consid- 
erable interest for a long period. Since the discovery of 
the anomalous low-temperature behavior of amorphous 
dielectrics |lj there were many attempts to explain this 
unusual phenomenon. And yet, in 1986 Freeman and An- 
derson Q reviewing both numerous empirical results and 
theoretical models made a rather drastic conclusion that 
the thermal conductivity of amorphous dielectrics " is not 
understood in any temperature range" . While at first 
glance this conclusion looks too pessimistic, it reflects an 
actual situation existed in glassy physics at that time. It 
has been found experimentally that the low-temperature 
anomaly in thermal conductivity, k, observed by Zeller 
and Pohl for some noncrystalline solids |IJ is universal. 
In particular, a T 2 dependence of k below IK together 
with the following characteristic plateau region have been 
found in the most of glassy materials (see, e.g., review 
as well as in some quasicrystalline alloys [|J. This find- 
ing allows to assume that a realistic theoretical model for 
description of the thermal conductivity should capture 
the very essence of these materials, their microstructure. 
Thus the conclusion of Freeman and Anderson actually 
suggests that this problem still remains to be solved. 

The most widely used model which explains glassy 
physics below IK has been proposed by Phillips || and 
Anderson, Halperin, and Varma In this tempera- 
ture region, thermal transport takes place via propagat- 
ing acoustic modes. Thus, k depends on the scattering 
of the phonons by the structure. Within the model H[|, 
the principal scatterers were proposed to be the tunnel- 
ing states (TS) in the glass. However, in spite of the 
success achieved in interpreting experimental data, some 
important questions were left unanswered. The main of 
them concern the nature of TS and the source of their 
universality (see, e.g., discussion in j||7|]). 

Since 1972 there have been many attempts to ex- 
plain the experimental data in terms of propagating 
thermal phonons scattered by TS |||§-|ll|] (in different 
modifications), density fluctuations (Rayleigh scatter- 
IjIQUQ, and the like. Some papers invokec 



ing) 



nonpropagating vibrational modes as well Sj. The most 
successful approach includes a combination of TS and the 
Rayleigh scattering JTG] . Notice that there were attempts 



to go beyond the TS approximation (see Ref. |H| and the 
references therein) as well as to do without TS at all. For 
example, an interesting phenomenological approach was 
discussed in Ref. |3j where a combination of the structure 
scattering together with a constant free path was shown 
to give a good fit to the thermal conductivity data in 
amorphous arsenic (a,- As). As a result, it was suggested 
that no tunneling state scattering is needed to explain 
the temperature dependence of k in this material. 

Of our special interest is the model proposed by J.C. 
Phillips |ll|, who suggested that glasses have cluster 
structures. It has long been known that there is a 
close connection between grain boundaries in materials 
with pronounced cluster structure and wedge disclination 
dipoles (WDD) J!?]]. Recently, this analogy has been suc- 
cessfully applied for description of the old experiments 
in polycrystals on the basis of the WDD-induced scat- 
tering [p"8| . As was mentioned in B, the scattering in 
polycrystals has a close parallel in amorphous solids. 

Another approach which takes into account topological 
defects in glasses came from geometrical and topological 
considerations. For example, Kleman (see review [fi"9"|) 
proposed that disordered covalent systems can be de- 
scribed in terms of distributions of defects, most prob- 
ably disclinations, in a crystal situated in a space of con- 
stant curvature. This idea was explored by Nelson |2C| ] 
in description of metal glasses. Similarly, Duffy and Riv- 
ier pH considered a continuous random network formed 
by 5 and/or 7 rings in the hexagonal lattice. In this 
case one can also expect the appearance of 5-7 disclina- 
tion pairs (dipoles of disclinations). Let us mention also 
papers where elastic dipoles have been introduced 
within the glassy model. 

In this Letter, we present a new look at the problem 
of thermal conductivity in dielectric glasses as well as 
at their microstructure. Actually, we reexamine the old 
ideas of the previous two paragraphs from another point 
of view. The model proposed here was inspired by a re- 
cent finding that a biaxial WDD is a rather unique scat- 
tcrer |l^] . Namely, it has been shown that the mean free 
path Id of thermal phonons scattered by biaxial WDD 
with nonskew axes of rotation behaves like Id ~ <? 
at small q while Id —> const with q increasing. As is 
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known, a similar behavior, but with l q ~ q for small q, 
has been found in old experiments on glasses J2^]. How- 
ever, the later experiments showed that in fact l q ~ q^ 1 
in agreement with the behavior of Id- The main aim 
of this paper is to show that the phonon scattering due 
to biaxial WDD, being combined with the Rayleigh scat- 
tering, can explain the thermal conductivity experiments 
in amorphous dielectrics over a wide temperature range. 
As an example, we compare the results obtained from 
the experimental data for vitreous silica (a-S^C^)- More- 
over, considering WDD as the basic structural elements 
of glasses one can explain the universal behavior seen in 
different amorphous dielectrics. 

To gain a better understanding, let us discuss briefly 
the generally accepted approaches as well as the main 
problems pointed out in H. The starting point for the 
calculations is the known kinetic formula 



C{uj,T)vl(uj,T)duj, 



(1) 



where C(uj,T)duj is the specific heat contributed by 
phonons within the frequency interval cko, v is an average 
phonon velocity, and 1(uj,T) is the mean free path of a 
phonon. As usual, at low temperatures we can assume a 
constant v and a density of phonon states quadratic in to 
with the Debye cutoff ujd- To reproduce the observed k 
over the full temperature range the best form of the total 
mean free path was taken to be [fill 



l(u; : T) = (l- 1 + l R 1 r 1 +lr, 



(2) 



where It is the mean free path due to scattering from 
TS's, Ir comes from the Rayleigh scattering, and l m in is 
constant which must be included to prevent I from be- 
coming unphysically small. It should be stressed that it 
is this required cut-off l m in that causes the main prob- 
lem. In fact, it has no natural explanation. A differ- 
ent situation arises in a- As where the experimental data 
were described by a combination of a constant mean free 
path resulting from phonon scattering by small holes in 
samples and the Rayleigh type scattering. In this case, 
the Rayleigh scattering becomes important in the high- 
frequency limit. Thus, one has to use the interpolation 
formula which describes the Rayleigh scattering over the 
complete frequency range. It reads H 



h 



(3) 



where D is a constant which has been considered as a 
fitting parameter, and fc^ is the Boltzman constant. This 
approach explains the thermal conductivity over the full 
temperature range without the TS scattering. However, 
this model is restricted to description of a,- As where k ~ 
T 3 has been observed below 0.5K instead of k ~ T 2 for 
the most of glasses. 

Let us return to the paper || where the following 
three characteristic temperature regions have been dis- 
tinguished: 



1. At T/Q < 0.01 the phonon scattering is such that 
l/X w 150 for any type of glass. Here O is the Debye 
temperature, and A is the wavelength of a phonon. The 
main question is why this relation is valid for such a var- 
ied spectrum of glasses? 

2. In the plateau region, 0.01 < T/Q < 0.1, the mean 
free path of the phonons rapidly decreases with increasing 
q. What is the scattering mechanism responsible for this 
frequency dependence? 

3. At T/Q > 0.1 it is not even clear what excitations 
are responsible for thermal transport. 

We will try to answer some of these questions below 
within our model. First of all, let us reproduce briefly 
the main results of the previous paper where the scatter- 
ing of thermal phonons by WDD was studied For 
simplicity, in what follows we consider only biaxial WDD 
with nonskew axes of rotation (below the term "biaxial 
WDD" will refer only to this type of biaxial WDD). 

As is known, dislocations and disclinations cause ad- 
ditional strain fields which lead to the scattering of 
phonons. A mean free path of phonons of frequency u> 
scattered by the potential associated with a static de- 
formation of a lattice caused by straight WDD can be 
calculated within the generally accepted approach (see., 
e.g. p3,55|). An effective perturbation energy due to 
the strain field caused by a single WDD is written as 
U(r) = TujJ^/SpEab, where fiu is the phonon energy with 
the wavevector q, to — qv, 7 is the Griinciscn constant, 
and Eab is the strain tensor due to WDD. For a biax- 
ial WDD, the strain matrix is found to have a complex 
form. To simplify the problem, we assume that incident 
phonons lie normally to disclination lines and choose the 
suitable geometry with disclination lines directed along 
the z-axis and the dipole's arm oriented along the x-axis. 
In this case, U(r) takes the form: 



B ( x + L) 2 +y 2 

U(x,y) = — In- 

2 (x — Ly + y z 



(4) 



where B = hqv^vil — 2cr)/(l — a), v is the Frank index, a 
is the Poisson constant, and 2L is the dipole separation. 

In accordance with Eq.(4), the problem reduces to the 
two-dimensional case with the phonon mean free path 
given by 



2tt 

C 1 =rii / (1 -cos 6)U(6)d6, 
'0 



(5) 



where 31(6) is an effective differential scattering radius, 
and iii is the areal density of defects. Within the Born 



approximation, K(#) is determined as 25 



&(0) = 



qS 2 



2-kTiv 2 



<q\U(r)\q>> 



(0) 



where all vectors are two-dimensional ones, p = q — q 1 ', 
S is a projected area, and the bar denotes an averaging 
procedure over a which defines an angle between p and 
the x-axis. In other words, it means the averaging over 
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randomly oriented dipoles in the xy plane. Evidently, the 
problem reduces to the estimation of the matrix element 
in Eq.(6) with the potential from Eq.(4). The result is 

AniB 

U (p, a) = 7. — sin(pi cos a). (7) 

After averaging of \U(p, a)\ 2 over a and following inte- 
gration in Eq.(5) with respect to 6 one obtains 

= %A 2 {vL) 2 n t q( K Jl{2qL) + j\ (2qL) - 

^ 1 -J (2qL)J 1 (2qLfj, (8) 

where A = irj(l — 2er)/(l — er), and J n (t) are the Bessel 
functions. Notice that the behavior of Id in Eq.(8) is 
actually governed by the only parameter 2L which char- 
acterizes the dipole separation. In the long wavelength 
limit one gets Id ~ q^ 1 while for A < L we obtain that 
Id — » const. It should be stressed that the appearance 
of the q-independent region distinguishes remarkably this 
scatterer from the others known. 

The total phonon mean free path is obtained by com- 
bining the inverse lengths 

ltot(u) = {lR 1 +l D 1 )- 1 - (9) 

with In and Id from Eq.(2) and Eq.(8), respectively. An 
important consequence of Eq.(9) is that we do not need 
any cut-off procedure in calculating k. 



the WDD-induced scattering, then there appears a re- 
gion with the two scattering processes involved, and fi- 
nally the Rayleigh scattering becomes dominant at high 
frequencies. This specific behavior results from the fact 
that both In and Id tend to constants when q increases. 
Obviously, any q dependence of scatterers would break 
this picture at large q. The main parameters which ac- 
tually define Itot are D from the Rayleigh scattering and 
2L from the WDD scattering. 
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FIG. 2. Phonon mean free path hot as a function of fre- 
quency at 2L — 20A (solid line), 2L — 40A (dashed line), and 
2L = 100 A (dotted line). The parameter set is the same as 
in Fig.l. 
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FIG. 1. Phonon mean free paths Id (dashed line), Ir 
(dotted line), and hot (solid line) as functions of frequency. 
The parameter set used is: rn — 5 x 10 11 cm -2 , 2L = 20A, 
v = 0.1, v = 4.1 X 10 5 cm/sec, A = 2.6, D = 1.0 cm _1 Jl" 4 , 
and l = 15A. 



As an example, Fig. 2 shows Itot for different dipole 
separations. It is interesting to note that the size of the 
transition region decreases remarkably with 2L increas- 
ing. Notice that this region alone is responsible for the 
plateau in the thermal conductivity. 
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Fig.l shows Id, Ir, and Itot for a characteristic set of 
the model parameters. Notice that the parameters for the 
Rayleigh scattering are chosen to be in agreement with 
those for &-SiC>2 |§ • One can clearly see that at low fre- 
quencies the dominant scattering mechanism comes from 



FIG. 3. Thermal conductivity vs temperature calculated 
according to Eq.(10) with l to t from Eq.(9) at lo = 10 A (solid 
line) and lo = 15A (dotted line). The parameter set is the 
same as in Fig.l, Q = 342K. Experimental data for &-Si02 
from Ref.[l] are indicated by triangles. 
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The last step is to calculate k in Eq.(l) with l tot . For 
this purpose, it is convenient to rewrite Eq.(l) in the 
dimensionless form 



2ir 2 h 3 v 2 



e/r 



x 4 e x {e x - 1) 2 l tot (x)dx, 



(10) 



where x — ftu/kBT, and the specific heat capacity is 
chosen in the standard Debye form. The result is shown 
in Fig. 3 where the experimental points for &-SiC>2 are 
marked by triangles. As is seen, there is a good agree- 
ment over a wide temperature range. For l = 15 A some 
deviations from the experimental data appear above AOK 
but the correct qualitative behavior retains even at high 
temperatures. The best fit for a-SiO^ is obtained for 
2L = 20A. Notice that this value agrees with a cluster 
size expected in vitreous silica |l6| ]. Let us stress once 
more that we have no need for poorly defined parame- 
ters like the cut-off l m in m Eq.(2) or an additional to In 
constant free path proposed in Q to fit the data below 
0.57^. 

Now we can try to answer the questions formulated in 
the beginning of this paper. In accordance with Fig.l, at 
low frequencies the dominant contribution comes from 
Id- One can easily estimate the value of Id/X in this 
region. Spreading out Eq.(8) at qL <C 1 one gets 



D 



1 



A 8irA 2 (vL) 2 m' 



(11) 



This is the constant which depends on the model pa- 
rameters. Thus, in accordance with our model, the only 
parameters which determine the behavior of glasses at 
low temperatures come from their microstructure which 
includes biaxial WDD as the basic elements. The main 
parameters are the density of WDD, their Frank index, 
and the dipole separation. In particular, for our choice 
of parameters, ^_d/A ~ 135. Let us emphasize that the 
experimental data for various glasses given in |lj] can be 
described by minor variations of these main parameters 
(together with D). Thus the relation ^tot/A ~ 150 is in- 
deed satisfied for all these dielectric glasses. Based on 
the results of our analysis, we can conclude that: 

(i) the WDD-induced scattering of phonons can explain 
the established relation l/X ~ 150 for various glasses; 

(ii) in the plateau region both mechanisms (WDD- 
induced and the Rayleigh scattering) should be ac- 
counted for to produce a necessary behavior; 

(iii) with increasing temperature, the Rayleigh scatter- 
ing becomes dominant. While there are some deviations 
from the experimental data at high temperatures, the 
general behavior looks quite correct. Thus nothing else 
but phonons transport the heat at T > 0.10, and no any 
other excitations are needed. 

To summarize, we have presented a new look at the 
problem of the thermal conductivity in amorphous di- 
electrics. The main suggestion is based on the intro- 
duction of biaxial WDD as the basic structural elements 
of glasses. The unique properties of the WDD-induced 



phonon scattering enabled us to describe the universal 
behavior of the thermal conductivity of amorphous di- 
electrics over a wide temperature range. Thus, the re- 
sults obtained support the cluster picture of amorphous 
dielectrics suggested by Phillips |H| and shed some light 
on the relationship between the microstructure of glasses 
and their transport properties. 
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